Competing ferromagnetism in high temperature copper oxide superconductors 
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The extreme variability of observables across the phase diagram of the cuprate high 
temperature superconductors has remained a profound mystery, with no convincing 
explanation of the superconducting dome. While much attention has been paid to the 
underdoped regime of the hole-doped cuprates because of its proximity to a complex 
Mott insulating phase, little attention has been paid to the overdoped regime. Ex- 
periments are beginning to reveal that the phenomenology of the overdoped regime 
is just as puzzling. For example, the electrons appear to form a Landau Fermi liquid, 
but this interpretation is problematic; any trace of Mott phenomena, as signified by 
incommensurate antiferromagnetic fluctuations, is absent, and the uniform spin sus- 
ceptibility shows a ferromagnetic upturn. Here we show and justify that many of these 
puzzles can be resolved if we assume that competing ferromagnetic fluctuations are 
simultaneously present with superconductivity, and the termination of the supercon- 
ducting dome in the overdoped regime marks a quantum critical point beyond which 
there should be a genuine ferromagnetic phase at zero temperature. We propose new 
experiments, and make new predictions, to test our theory and suggest that effort 
must be mounted to elucidate the nature of the overdoped regime, if the problem of 
high temperature superconductivity is to be solved. Our approach places competing 
order as the root of the complexity of the cuprate phase diagram. 



The superconducting dome (See Fig. [T]), that is, the 
shape of the superconducting transition temperature Tc 
as a function of doping (added charge carriers), x, is 
a clue that the high-Tc superconductors are uncon- 
ventional. Conventional superconductors, explained so 
beautifully by Bardecn, Cooper, and Schrieffer (BCS) 
have a unique ground state that is not naturally sepa- 
rated by any non-superconducting states. The electron- 
phonon mechanism leads to superconductivity for arbi- 
trarily weak attraction between electrons. To destroy a 
superconducting state requires a magnetic field or strong 
material disorder. In the absence of disorder or mag- 
netic field, it is difficult to explain the sharp cutoffs at 
xi and X2 within the BCS theory. For high-Tc supercon- 
ductors, there is considerable evidence that competing 
order parameters are the underlying reason. Thus, xi 
and X2 signify quantum phase transitions, most likely 
quantum critical points (QCP). [1] Understanding high- 
Tc superconductors therefore require an understanding 
of possible competing orders. [1,|3| The QCP at Xi is ex- 
tensively studied, [4, 4 S S [1| but little is known about 
X2- Recent work has also emphasized the importance of 
the maximum of the uniform susceptibility in defining the 
pseudogap line T* in Fig. [T] that ends at another QCP at 
Xc- Q Here, we attempt, instead, at gaining insight from 
the possible existence of a QCP at X2- 

Complex materials — cuprates, heavy fermions, 
organics — are rife with spontaneously broken sym- 
metries in their ground states. This is unsurprising 
from a theoretical perspective. Because two distinct 
macroscopic states cannot be analytic continuations of 
each other, they cannot be literally said to arise from 
each other. Macroscopic properties are best described 



by a low-energy effective Hamiltonian. In the process 
of coarse-graining, one generates effective Hamiltonians 
with a multitude of coupling constants that signify 
the intermediate processes removed from the original 
Hamiltonian. Tuning these coupling constants can easily 
lead to a variety of ordered states that compete with 
each other. Thus, it is a truism that any symmetry that 
can be broken, must be broken. Competing phases can 
therefore explain the rich phase diagram of the high 
temperature superconductors. 

Nonetheless, it is frequently assumed that the super- 
conductivity in the overdoped cuprates (See Fig. [1]) 
grows out of a Landau Fermi liquid. Experimental sup- 
port has included the observation of angular magnetore- 
sistance oscillations in Tl2Ba2Cu06+5 (T12201) ^ and 
the reported verification of the Wiedemann-Franz law 
(WF) in both T12201 and La2-xSrxCu04 (LSCO), 0, 
lll| materials that can be overdoped to the edge of the 
superconducting dome. The Fermi liquid interpretation 
has not been without difficulties, however. 1^ The ver- 
ification of WF in Ref. can be questioned, as the 
relation was found to coexist with a large linear resistiv- 
ity term extending down to zero temperature. Given a 
number of difficulties, this is not clear even for Ref. 
An example is the failure of the Kadowaki- Woods (KW) 
relation. [l3| The resistivity, p , of a Fermi liquid, no 
matter how strongly interacting it may be, is given by 
P = Po + AT'^, T ^ 0, where po is its residual value 
and T is the temperature; similarly, its specific heat is 
given by C = 'yoT, T ^ 0. The relation between A and 
7o is known as the KW relation, and it is easy to see 
that the KW ratio, A/'^q , should be independent of the 
effective mass to*, that is, the band mass to renormal- 
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FIG. 1: A sketch of a simplified phase diagram of the cuprates 
as a function of temperature, T, and carrier concentration 
(doped holes), x. The superconducting dome rises at about 
x\ « 5% and ends at about X2 ~ 30%; xi and X2 are two pro- 
posed QCPs terminating the d-wave superconductor (DSC). 
The regime proximate to a::i is termed underdoped and that 
close to X2 is designated overdoped. The enigmatic pseudogap 
state is marked by the line T* . A plausible candidate for this 
state at low temperatures is the d-density wave (DDW) i. 
The doping Xc marks another QCP, and roughly the region 
between Xc and X2 is the overdoped regime. The underdoped 
regime is riddled with other non-superconducting competing 
phases, such as stripes, checkerboard charge order, spin glass, 
and the three-dimensional antifrromagnet (3D AF). Here X2 
marks the formation of a ferromagnetic (F) phase. The quan- 
tum critical regime emanating from X2 is labelled QCR. The 
arrows (a), (b), and (c) are possible experimental trajectories. 
It is only for (c) that a static ferromagnetic order parameter 
can develop 



ized by interactions. The striking enhancement ^ of 
the KW ratio in heavily overdoped LSCO and T12201 
over the "universal" heavy fermion (to* 3> m) value, but 
with TO* « TO , is surprising. The anomalous character of 
the overdoped Lai.7Sro.3Cu04 is also clear [11] from the 
extremely narrow range of temperature for which the lin- 
earity of specific heat holds, while the quadratic temper- 
ature dependence of the resistivity holds over a surpris- 
ingly broad temperature range up to 50 K, beyond which 
it behaves as T^-^ . There is little doubt that electron- 
electron interaction cannot be subsumed by Fermi liquid 
renormalizations, and the measured enhancement of the 
KW ratio is not particularly meaningful. Clearly, further 
experiments are necessary. 

The centerpiece of the anomalies is the uniform mag- 
netic susceptibility, which does not remotely resemble 
that of a Fermi liquid. It is known that in La2-xMxCu04 



(M=Sr, Ba), the susceptibility exhibits a dramatic up- 
turn at temperatures just above the transition temper- 
ature Tc for X > 0.2, while for lower doping it behaves 
as it would if the ground state were a singlet. Similarly, 
overdoped T12201 shows a striking upturn in suscepti- 
bility; see Fig. [51 In both cases, the upturn becomes 



14 



more pronounced with overdoping. 
Despite motivated effort, it has not been possible to at- 
tribute this unusual behaviour to Curie susceptibility of 
a magnetic impurity phase. An explanation in terms of 
localized moments, as suggested by Nakano et al. [20| . 
is also problematic, because the effective Curie constant 
would have to increase with hole doping, when the sys- 
tem becomes more metallic. This seems very unlikely. 
Moreover, a careful analysis shows that 1/4 of Sr ions 
in La2-xSrxCu04 corresponding to holes that exceed 
X = 0.22 have to create paramagnetic moments, which 
is equally unlikely. [l3| Neutron scattering measurements 
around q = (tt, tt) in the same sample show anomalous 
enhancement of dynamic antiferromagnetic correlations 
due to Zn substitution, while Sr overdoping results in a 
smooth decrease of the imaginary part of the dynamic 
susceptibility. Given that Zn substitution is known pro- 
duce paramagnetic moments, the different behavior of Sr 
overdoping argues against the formation of paramagnetic 
moments. 

One of the hallmarks of special correlation effects that 
give rise to Mott insulators (as opposed to band insu- 
lators) is antiferromagnetism, which takes on an incom- 
mensurate character with doping. Recent neutron scat- 
tering experiments in LSCO indicate, however, that the 
intensity of the finite energy incommensurate antiferro- 
magnetic peaks diminishes with increased doping, van- 
ishing as the material is tuned to the edge of the su- 
perconducting dome at x = 0.3. [21[ This is a striking 
testament to the disappearance of these correlations in 
the overdoped regime. 

Here we suggest that ferromagnetic fluctuations com- 
pete with d-wave superconductivity in the overdoped 
regime. Ferromagnetism can be a consequence of strong 
electronic correlations and low particle density (large hole 
density) when Fermi surface nesting is remote or the 
system is not close to a Mott insulating phase, [l^, [2^ 
Anderson [23| has also noted the possible role of ferro- 
magnetic correlations in the overdo ped regime from the 
perspective of Nagaoka mechanism 25|, |26| approaching 
it from the underdoped side. We consider, instead, the 
itinerant picture and use it not to justify the pseudogap 
crossover, but the rich phenomenology of the overdoped 
cuprates. It would be interesting to study the intersec- 
tion of the two ideas. We suppose that X2 (See Fig. [T]) 
signifies either a QCP or a first order quantum phase 
transition between a d-wave superconductor and a fer- 
romagnetic ground state. A first order transition can 
become continuous in the presence of disorder — for or- 
der parameters of continuous symmetry, this is expected 
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FIG. 2: Total magnetic susceptibility in a field of IT of T12201 
in the overdoped regime with Tc tuned to zero (data excerpted 
from Ref. [13]). The data can be fit to the form x = A + B/T'< 
globally over the entire range. The core diamagnetic suscepti- 
bility, — 1.93 X 10~*emu/mol, is independent of temperature. 
The range between 30 K and 100 K can be equally well fit 
(within the accuracy of the data) with 7 between 0.95 — 1.25. 
The inset shows a fit with 7 = 5/4, a prediction of the quan- 
tum critical scaling given in the text. Clearly quantum criti- 
cal scaling should become more precise at asymptotically low 
temperatures for which data are not available. See Ref. [13 ] 
for other hole concentrations. Similar data are also available 
for LSCO and LBCO (See Refs. [H,!!!). 



to occur in dimensions d < 4, |27|, |28|, m, [3^, im which 
can be straightforwardly extended to a quantum phase 
transition between two ordered states. Thus, ferromag- 
netic fluctuations in the superconducting state will act as 
pair breakers and explain the demise of superconductiv- 
ity in overdoped cuprates. Of course, pair breaking will 
occur even if the transition is weakly first order — large 
ferromagnetic domains will still appear in the supercon- 
ducting phase. 

The possibility of the coexistence of two phases would 
have been plausible if the superconductivity were of the 
triplet variety. In principle ferromagnetic fluctuations 
can give rise to spin-triplet p-wave supgrconductivity, but 
this appears to be far too delicate. [3^ In any case, at 
least approximately up to 22% doping, phase sensitive 
experiments seem to indicate robust d-wave pairing in 
the cuprates fssi, which is also generally believed to be 
true within the entire superconducting dome. 

It was argued elsewhere that the criticality at xi is 
that of a quantum three-dimensional Xy-model, [l| the 
fundamental variable being the phase of the supercon- 
ducting order parameter. This alone leads to the robust 
result Tc oc ps (0)1/2^ where psiQ) is the T = super- 
fluid density. [6| Recent experiments 34] are consistent 



with this quantum critical result asymptotically close to 
xi , showing, a strong deviation from the widely assumed 
Uemura plot, Tc oc Ps(0). [s^ 

What would be the corresponding result at X2^ if it 
is a QCP as well? For ] a; — 0:2] ^ 1, pair breaking will 
lead to Tc cx \x — X2\^ , where is an exponent to be 
determined. Because of spin diffusion, the dynamic ex- 
ponent, 2i for a dirty metallic ferromagnet is given by 
2 = 4, [Sy] dissipation being relevant on both sides of 
the transition; on the superconducting side it arises from 
the gaplessness. From an interesting experiment 37| on a 
classical system, it appears that a disorder driven contin- 
uous transition can be of emergent 0* -type (Ising) , [sl] 
while the actual phases on either sides of the transition 
are determined by the operators that are irrelevant at the 
critical point. We surmise that the same is true for the 
quantum problem, except for an added dissipative term 
in the effective action. 3^, 4^ Since for d — {d+z) > 4, 
Ps(0) ]a: — 3:2], the mean field result. Due to the danger- 
ous irrelevant operator, u , controlling the interaction of 
the bosonic modes, a one-loop calculation (higher loops 
cannot change the leading scaling result) leads to the 
correlation len gth ^(T, x = 0:2) ~ l/T^^+^^Z^^, consis- 
tent with Ref. [40|. But as the phase boundary is given 
by ^(T = Tc,x = X2) ^ CiT = 0,x), the exponent 9 
is given by 6* = z/(z -I- 1), using the mean field result 
^(T = 0,a;) ~ \x — X2\^^^^- If we substitute z = 4, we get 
9 — 4/5 and Tc oc ps(0)^/^, in contrast to the underdoped 
regime, providing an explanation of the Boomerang effect 
of Tc vs. Ps{0) discussed earlier, [41[ where pair breaking 
of unspecified origin was suggested. This exponent can 
be tested in future experiments. (The exponent in the 
overdoped regime derived in Ref. [6| would be 2, which 
is incorrect.) 

Recently an empirical scaling analysis [i^] shows that 
Ps(0) oc adciT « Tc)Tc holds reasonably well for a broad 
range of doping, both for in-plane and out-of-plane com- 
ponents, where the quantity (Tdc(T w Tc) is the conduc- 
tivity at r « Tc. A similar scaling relation was derived 
earlier for the out-of-plane components. Ja^I Even if this 
relation were to hold for the entire doping range, this 
would not invalidate the Boomerang effect, which is the 
statement that the T — exponents of Ps(0), as a func- 
tion of X, are different at the QCPs xi and X2- This 
would merely mean that these exponents are the sum 
of the exponents of the conductivity and the transition 
temperature as a function of x at the respective QCPs. 

Ferromagnetic fluctuations result in the surprising up- 
turn of the susceptibility. Previous fits to the data have 
assumed a 1/T-dependence. For a first order transition, 
this may be an adequate description of the contribution 
from finite-size ferromagnetic domains. But if a QCP ex- 
ists at X2, the susceptibility, x, will exhibit quantum crit- 
ical scaling; for x tuned to X2, x ^ S,{T,x — X2Y ^ T^'' , 
as T ^ 0, where 7 = 5/4 (z = 4). This is consistent with 
the data shown in Fig. [2] for Tl 2201, which, to a good 
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approximation, is tuned to X2. Note that x approaches 
the Curie form as z — > cx) and — > 1. Importantly, for 
the range of doping explored in experiments thus far, the 
data cannot, even in principle, reflect a static ferromag- 
netic order; for this, the doping must be larger than X2 
(See Fig. [1]). From a one-loop theory (similar to the sus- 
ceptibility calculation) we predict that for x = X2 the 
magnetization, M, and the applied field, h, must obey 
(also approximately true in the entire quantum critical 
region) : 



M = {h/uY^'f{uh\') 



(1) 



where the universal scaling function f{y) — > 1 as y — > oo, 
and f{y) y^^^ as y ^ 0. Here d = 3 and {d + z) > 4. 
The crossover from the finite temperature result to the 
zero temperature result is given by ~ {h/uY^^. Since 
presently there are no known samples for a; > 2:2 , we can 
only verify this quantum critical scaling by either suscep- 
tibility measurements down to low temperatures, or by 
precise field dependence of Knight shift, a measurement 
that to our knowledge has not been carried out. 

It is useful to sketch this calculation in a bit more 
detail. Consider a iV-component (j)* theory with dissi- 
pation corresponding z = 4 such that {d + z) > 4, that 
is, above its upper critical dimension. We integrate out 
the Matsubara frequencies 7^ to construct an ef- 
fective d-dimensional theory, [ll, [s^ The leading scaling 
result can be obtained from a one-loop calculation; higher 
order calculations can not change the leading behavior. 
The spatially independent part of the quadratic term, R, 
and the coefficient of the quartic term, U, are 
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U = u + Oiu'^), 



x{q,iLUn), (2) 
(3) 



where x~"'^(q, icj„) = q'^ + \uJn\/Dq^ + r, and r oc (2:2 — 
x) is the original quadratic term in the quantum {d + 
z)-dimensional problem. Here D is the spin diffusion 
constant of an itinerant metallic ferromagnet. There is 
of course no renormalization of the term containing the 
static magnetic field h. Tuned to the T = quantum 
critical point at r = 0, the equation for the magnetization 
becomes, 



C'^{r = 0, T)M -h+ = 



(4) 



where 



r'(r = o,r) = 



N- 



dq 



^11 / 7^^o(q'*^»)' 
(5) 

and Xo^i^i^^'-^n) = + luJul/Dq"^- Note that the sum 
now includes n = 0. The solution of the cubic equation 
in Eq. |3]now gives the universal scaling function in Eq. [T] 




FIG. 3: The local density of states (LDOS) (arb. units) as a 
function of frequency to, N{r,u}), in a d-wave superconductor 
(DSC) with a central ferromagnetic domain of size 5x5 em- 
bedded in a sample of size 41 x 41 (in terms of lattice units). 
The curves labelled (m, n) correspond to the location where 
the density of states is calculated from Bogoliubov-de Gennes 
theory. The density of states at (0, 0) consists of an extraor- 
dinarily sharp peak with a small gap. As we move out of 
the centre, these peaks decrease in intensity, get wider, and 
eventually the LDOS corresponding to the surrounding DSC 
take over. This effect is demonstrated by calculating LDOS 
at (2.0) and (41,41). We see sharp peaks in the central re- 
gion even though the depression of the order parameter due to 
ferromagnetism is only 5%. Note also the change in the low 
energy spectra as we move out of the central region, which 
becomes linear, as it should for a pure DSC. The parameters 
used in this calculation are listed in the methods section. 



As to excitations in the superconducting state, quan- 
tum critical fluctuations will primarily scatter the nodal 
quasiparticles as x — > X2, because their energy ap- 
proaches zero; in degenerate perturbation theory, modes 
of equal energies couple most strongly. Thus, we expect 
nodal quasiparticles to broaden with overdoping while 
the antinodal quasiparticles, whose energy is finite, will 
sharpen, because they will not couple strongly to the 
critical fluctuations. This surprising effect may have 
already been observed in recent angle resolved photoe- 
mission measurements. 



44l | This is in part due to the 



remarkable effort in recent years in the preparation of 
high quality single crystals of overdoped T12201. [45*1 At 
present two models involving impurity scattering may 
be relevant, [i^, [13] although these predict a isotropic 
quasiparticle scattering rate above Tc and a temperature 
dependence set by the scale of the energy gap. More 
experiments are necessary to discriminate between dif- 
ferent ideas, but such substantial impurity scattering in 
the high quality single crystals seems unlikely. As to in- 



5 



elastic mechanism, another possibility is the proximity of 
a competing dx2_y2+idxy superconducting phase, similar 
to our suggestion of a competing ferromagnetic state, [i^ 

The STM (scanning tunnelling microscopy) spectra 
should be sensitive to the inhomogeneity due to ferro- 
magnetic domains. We used a self-consistent Bogoliubov- 
de Gennes theory to calculate the local density of states 
(LDOS) for a d-wave superconductor modified to include 
an effective Zeeman field. The calculational method has 
been described in detail in the past. [4§| Here we simply 
note the modifications necessary for the present problem. 
The Hamiltonian should contain an added ferromagnetic 
term Hp ^ ^JFj2{ij) ' where the sum is over 
nearest neighbour pairs of spin operators St . We approx- 
imate this term in mean field theory as a Zeeman field 
—h Sf ■ While h is treated as a fixed quantity, all other 
parameters are determined self-consistently. The remain- 
ing Hamiltonian is a t — i' — J- model, where the units are 
such that t = I, t' = -0.27. J = 1.27. These are ef- 
fective parameters including mean field renormalizations 
(Gutzwiller-constraint) due to hole doping oi 5 = 25%. 
These parameters result in a uniform DSC with a d-wave 
gap A/j = (Ao/2){cosfc:r — cosfc^), where Aq = 0.44. 
For Fig. [HI we chose /i/Ao = 0.4; the reduction of the 
self-consistent pairing amplitude at the centre of the fer- 
romagnetic patch with respect to Aq was about 5%. The 
size of the unit cell was 41 x 41 replicated in the repeated 
zone scheme 14 x 14 times. The central ferromagnetic 
patch was of size 5x5. We have repeated the calcu- 
lation with a variety of parameters. The results shown 
in Fig. [3] are representative. The ferromagnetic exchange 
constant, Jp, corresponding to the assumed Zeeman field 
is much smaller than the ferromagnetic Jq that one needs 
to explain the phenomena in the overdoped regime in 
terms of magnetic impurity effects. The field, which is 
zero everywhere except for a small patch of sites at the 
centre of the lattice, represents a pair-breaking ferromag- 
netic domain. As shown in Fig. [3l the results of this cal- 
culation are consistent with the recent work where the 
gap in the central patch was arbitrarily set to zero. (50l | 
We see the same effect even though the amplitude is re- 
duced by only 5%, provided self-consistency is imposed. 

We emphasize, however, that the gap inhomogeneity 
in the underdoped regime [IH is entirely different. [4§| 
There the large gap regions reflect a pseudogap due to a 
competing orde — in our opinion, the d-density wave Q — 
while the narrow gap regions reflect the d-wave super- 
conductor. This change in character of the elementary 
excitations should take place around the middle of the 
dome (See Fig. [H). 

If ferromagnctism is truly a competing phase in the 
overdoped regime, the vortices would be easier to nucle- 
ate and the physics of the vortex core could be dramati- 
cally different, which could be probed in STM and other 
measurements. The most exciting test of our proposal 
will be a direct observation of the ferromagnetic phase, 



which would require the fabrication of samples with hole 
concentrations greater than X2- Since there can be no 
ferromagnetic order at flnite temperatures in two dimen- 
sions, the Curie temperature will be set by the interlayer 
coupling and is likely to be very small, especially for a 
disordered ferromagnet with direct exchange or dipolar 
coupling. Nonetheless, we hope that our ideas are suffi- 
ciently tantalizing to spark the interest of the experimen- 
talists. 
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